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ABSTRACT 

2 ■ Observations of microlensing transit events can be used to measure the limb 

O ■ darkening of the lensed star. We discuss the advantages and drawbacks of several 

' microlensing light curve inversion methods. The method of choice in this work is 

d • inversion by means of decomposition of the stellar surface brightness profile. We 

^ ■ construct an ideal basis by principal component analysis of brightness profiles 

Q> ' obtained from model atmosphere calculations. Limb darkening approximations 

^ ■ using such a basis are superior to those using standard power-law limb darkening 

^ ■ laws. We perform a full analysis of simulated single-lens microlensing transit 

^ ' events including a detailed error analysis of the method. In realistic events with 

^ ■ a low impact parameter the brightness profile of the source can be recovered with 

ir^ • a relative accuracy of 2% from the center of the source disk to 0.9 of the disk 

^ ■ radius. We show that in the particular case of the observed MACHO Alert 95-30 

O ■ event the intrinsic complex variability of the lensed red giant hinders efforts to 

^ ■ recover its surface features. 

2 ! Subject headings: gravitational lensing — methods: data analysis — stars: at- 

'. mospheres 



1. INTRODUCTION 

Our current knowledge of the intrinsic physical properties of stars is largely based on 
stellar atmosphere models (e.g., Gray 1992). Ever more advanced models are computed for 
a broad range of stellar types, taking into consideration a range of physical and chemical 
processes expected to occur in the star's atmosphere. 

However, for the vast majority of observable non-variable stars the predictions of these 
models can be confronted with a single observable quantity - the integrated full-disk spec- 
trum. It follows that any two physically different atmosphere models which predict similar 
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integrated spectra cannot be observationally distinguished. More specifically, any obser- 
vationally "verified" model may have a center-to-limb variation of the spectrum (or limb 
darkening, in terms of broad-band photometry) very different from the actual star and thus 
be unphysical, yielding wrong physical parameters of the star. This inherent problem can be 
overcome by testing atmosphere models in cases or events when the limb darkening of the 
observed star can be measured. 

Until recently the only stars other than the Sun providing the opportunity to directly 
study their surfaces were the components of eclipsing binaries. Measuring their limb dark- 
ening profiles from eclipse light curves proved to be no easy task, due to the degeneracy 
with other fitted parameters (Grygar, Cooper, & Jurkevich 1972). Even obtaining useful 
constraints on the value of the linear limb darkening parameter requires high quality hght 
curves (Popper 1984; Popper et al. 1985). 

Extra-solar planets occulting their parent stars are analogous to eclipsing binaries with 
small secondaries. High precision observations of these scarce systems can also be used for 
measuring the parent star limb darkening. This was recently demonstrated in the case of 
HD 209458, the first observed extra-solar occulting system (Jha et al. 2000; Deeg, Garrido, 
& Claret 2001; Brown et al. 2001). 

The surfaces of nearby giant stars can be resolved by interferometry. Developments 
in stellar interferometry have led to an advance from measuring stellar radii (Mozurkewich 
et al. 1991) to detecting the effect of limb darkening (Quirrenbach et al. 1996) to limb 
darkening measurement for a. Ori (Burns et al. 1997)^. Hajian et al. (1998) points out the 
main difficulty of such measurements - signatures of differences between limb darkening 
profiles become detectable only beyond the first zero of the visibility function. At these high 
spatial frequencies the fringe visibilities are low, and thus difficult to measure with sufficient 
precision. Moreover, as shown by Jacob et al. (2000), simultaneous fitting of the stellar 
radius and limb darkening is complicated by degeneracy. These and other studies indicate 
that while the use of interferometry for constraining the limb darkening of the source is 
promising, the accuracy will not reach the level achievable for eclipsing binaries in the near 
future. 

In this work we concentrate on another method - measuring limb darkening from ob- 
servations of source-transit microlensing events. Galactic gravitational microlensing occurs 
when a foreground massive object such as a dim star (hereafter "lens" ) approaches the line 
of sight to a background star ("source"). Due to gravitational deflection of hght by the 
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lens, the flux from the source is temporarily amplified as the lens passes in the foreground 
(Paczyhski 1996). Currently operating microlensing surveys have already detected over 500 
such events toward the rich stellar fields of the Galactic bulge and the Magellanic Clouds. In 
case the lens directly transits the source star, the light curve of the event will be affected by 
its limb darkening profile (Witt 1995; Valls-Gabaud 1995; Loeb & Sassclov 1995). Observa- 
tions of one such transit event, MACHO Alert 95-30, were reported by Alcock et al. (1997). 
By the time of writing, limb darkening parameters have been measured also in four binary- 
lens caustic-crossing events: MACHO 97-BLG-28 (Albrow et al. 1999a), MACHO 98-SMC-l 
(Afonso et al. 2000) ^ MACHO 97-BLG-41 (Albrow et al. 2000) and OGLE-1999-BUL-23 
(Albrow et al. 2001). ^ 

Microlensing thus presents a unique new possibility for studying the surfaces of distant 
stars (Valls-Gabaud 1995; Sasselov 1996). As discussed in §3, several approaches have been 
taken by other authors to extract the limb darkening profile from microlensing light curves, 
with encouraging results. Here we present another inversion strategy providing the most 
accurate and stable results so far. We achieve this by decomposing the profile using a basis 
constructed by principal component analysis of a set of actual model atmosphere profiles. 

In the following presentation, §2 contains the basic description of single point-mass 
lens microlensing of a stellar source. In §3 we discuss several methods of inverting the 
microlensing light curve to obtain the surface brightness profile of the source. Specifically, in 
§3.3 we provide the basic equations for inversion by brightness profile decomposition. Using 
model atmosphere calculations and principal component analysis of the model profiles, we 
determine an ideal basis for brightness decomposition in §4. The main results are presented 
in §5, in which we analyze simulated light curves, simultaneously recovering the lensing 
parameters and the brightness profile of the source. We perform a detailed error analysis 
and compute the inversion accuracy for a source with unknown limb darkening. We compare 
the results with those obtained by fitting the linear limb darkening law, applied in all five 
published events (an additional square-root coefficient was fitted in MACHO 97-BLG-28). 
The analysis of the observed MACHO Alert 95-30 event light curves is described in §6. Limb 
darkening measurement prospects are discussed in §7, followed by a summary of the main 
results in §8. Two appendices are added, the first describing broadband chromaticity of 
microlensing light curves, the second discussing the model limb darkening profiles and their 
fits. 



^This particular event yielded the first and thus far unique case of limb darkening measurement for an A 
star from a metal-poor population. 

^Preliminary limb darkening measurement was reported in a fifth event, EROS BLG-2000-5 (An et al. 
2002), during the review of this paper. 
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2. GRAVITATIONAL MICROLENSING LIGHT CURVES 

Throughout this work we describe the source and lensing geometry in terms of angular 
displacements in the plane of the sky, using the angular radius of the lensed star as a 
(projected) distance unit in this plane. ^ We concentrate here on the case of a spotless source 
star, with a symmetric surface brightness profile B{r). Microlensing of a spotted source is 
discussed elsewhere (Heyrovsky & Sasselov 2000). 

A point-mass lens located at a projected distance / from the center of the disk of a 
background star increases the observed flux from this source star to a value 



1 

F{1) = J B{r)A{l,r)rdr 



where r is the distance from the center of the star, and the surface brightness is measured 
in appropriate units. The angle-integrated amplification A{l,r) is defined by 



2tt 



A{l,r)^ J Ao{y/P - 2r I cos (j) + r^)d(f) , (2) 



where (p is the (source-centered) polar angle between a point on the stellar disk at a distance 
r and the lens at a distance I. The lensing amplification factor Aq{x) for a point source at 
a separation x from the lens is 

here e is the angular Einstein radius of the lens (Paczyhski 1996) in units of lensed star 
radius. The integral in equation (2) can be directly solved in terms of elliptic integrals as 
follows: 
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^Using the angular Einstein radius as a unit (traditional microlensing notation) is not convenient in 
this work concentrating on the source star structure. The conversions between our parameters and the 
parameters in the standardized notation proposed by Gould (2000) are I = uOe/O*, e = QeI^*-, V = uqOe/O* 
and i* = tE0*/9E- 
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We use standard notation for the complete elliptic integrals of the first (K) and third (11) 
kind. In the latter we follow the sign convention used by Byrd & Friedman (1971)^, 
U{a^ k) = - a2 sin2 e)-\l - k"" sin^ ^)-V2 dO . 

Notice that A{1, r) is symmetric in its arguments I, r. Also note that >A(0, r) = 27r74o(r), 
as well as yi(Z,0) = 2t^Aq{1). When the lens is far from the source {I » l,e), the angle- 
integrated amplification drops to 2tt. In the limit of a strong lens (e ^ 1), the angle- 
integrated amplification at the light curve peak (e ^ I) in equation (1) can be approximated 

by 



A{l,r) 




an expression linearly dependent on the Einstein radius e. In this work we use the more 
general expression in equation (4). Figure 1 contains sample plots of the angle-integrated 
amplification as a function of radial distance r for Einstein radii e — 1 and 10 (the latter 
scaled down by a factor of 10 for comparison), and several lens positions. At the position of 
the lens (r = I) there is a logarithmic divergence, so that lor r — I -\- 5 

where \5\ <^ 1. Due to this divergence at the lens position (for / < 1) in the integrand of 
equation (1), the lens efficiently scans the surface brightness profile of the source star as it 
moves across the stellar disk. The general microlensing effect on the spectrum of model stellar 
atmospheres was demonstrated by Heyrovsky, Sasselov, & Loeb (2000, hereafter HSL). We 
provide a brief discussion of broadband light curve shapes in Appendix A of this work. 



3. LIGHT CURVE INVERSION METHODS 

Photometric observations of microlensing events provide us with sets of observed fluxes 
Fi measured at different times ti with estimated errors cXi {i = 1, . . . n). Primary analysis of 
a transit event yields the main (geometric) lensing parameters: the Einstein radius e, impact 
parameter p, time of closest approach to and angular source radius crossing time i*. The 
latter three parameters help us mainly to convert the observation times ti to lens positions 
li , as shown further in equation (29). In the two following sections (through §4), we assume 
zero blended flux (contribution of unresolved stars) and consider the geometric parameters 
to be fixed at predetermined values, in order to facilitate the inversion analysis. In §5 and §6 



Press et al. (1992) use an opposite sign for the first parameter. 
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we proceed to the full event analysis including the effect of blending, and recover the lensing 
and limb darkening parameters simultaneously. 

The aim of this work is to extract the brightness profile B{r) from the observed fluxes 
Fj. The mathematical problem is to invert equation (1), classified as a Predholm integral 
equation of the first kind. This problem has no universal method of solution. The task is 
more difficult in this particular case, because of the required high level of accuracy of the 
result. The aim is not just to detect limb darkening, but to actually measure it. Based on 
solar observations and model atmosphere calculations, broadband limb darkening profiles in 
the optical range vary only to a fairly limited degree (see Appendix B), hence the required 
precision. 

There is a number of different numerical approaches to inverting equation (1) that can be 
tested. Craig & Brown (1986) provide a good overview of strategies for solving astrophysical 
inverse problems. In the following two subsections we discuss two methods applied previously 
by other authors to the microlensing case. The third subsection introduces the method used 
in this work. 



3.1. Inversion by Discretization 

A straightforward method is to invert equation (1) by discretization, approximating the 
integral by a sum. For a light curve with n points we can select ud < n radial points Vj 
(j = !,...«£)) on the source disk and write 

no 

5(r)yi(/„r)rrfr = ^D,,5(r,) (7) 



for each lens position (i = 1, . . . n). Here the matrix element ©jj is some approximation of 
J A{li, r) rdr over the vicinity of point rj. 

In order to obtain brightness point estimates B{rj) providing a good fit to the light 
curve data, we now minimize 

n 

x'-J2^-'[Fi-F{k)r . (8) 

i=l 

Rescaling the matrix to ©ij = ^ij/ui and the measured fluxes to Fi = Fi/ai, we can write 
the x^-minimizing solution explicitly: 

n no 

^(^.) = EE(^^^tf ^^^^^ ' (9) 

i=l k=l 
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here denotes the transpose of matrix D. The errors due to hght curve noise for each of 
the obtained brightness points are given by the standard deviations 



In comparison with the other methods discussed subsequently, this direct approach has the 
advantage of not making exphcit assumptions about the shape of the brightness profile. 

Gaudi & Gould (1999) presented an error analysis of this method®, demonstrating that 
the surface brightness profile can be recovered only for the outer radial part of the stellar 
disk directly transited by the lens. The profile of the inner part, at radii smaller than the 
impact parameter of the lensing event, cannot be accurately reconstructed. 

The general failure of microlensing light curve inversion in this region can be explained 
as follows. For any light curve point k > 0, the matrix element Djj — > for inner disk points 
rj — > by definition. Hence, in the case of a non-zero impact parameter the associated 
value B{rj) has negligible weight in all the sums given by equation (7), and thus is weakly 
constrained by the data. Alternatively, this failure can be understood by noting the low 
microlensing sensitivity to stellar spots located off the projected lens trajectory (e.g., a 
circular spot at the disk center; see Heyrovsky & Sasselov 2000). 

Returning specifically to inversion by discretization, some of its intrinsic problems are 
directly apparent from equation (7). For transit points (Zj < 1) the integrand is divergent, as 
demonstrated by equation (6). Obtaining an adequate approximation of such an integral by 
a sum requires a fairly high number of terms, i.e., an even higher number of transit light curve 
points. Even in the case of negligible light curve noise, simulations show that discretization 
with a limited number of points introduces errors that make the method unstable and the 
inversion unreliable. 

As a further point, straightforward discretization assuming constant surface brightness 
within annuli of the source disk can introduce additional errors. Unless there is a sufficient 
number of radial points defining the annuli, this assumption cannot be made particularly 
near the limb of the source, where the brightness drops rapidly. 

The stability of this method can be improved by introducing additional constraints on 
the solution points B{rj). For example, one might require positiveness of the solution, its 
monotonic decrease or concavity. Typically these constrain the solution points by a set of 
linear inequalities. The solution is then found numerically by minimizing given by equa- 



^Their paper covers the case of binary lensing as well, in the linear caustic crossing regime. Here we 
discuss their single lens results. 
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tion (8), while satisfying the imposed inequahties. Bogdanov & Cherepashchuk (1996) took 
this approach in the case of non-transit events^, seeking a non-negative concave decreasing 
solution. Unfortunately, they made no estimate of the errors or stability of their results. 
In a following paper (Bogdanov & Cherepashchuk 1999), the authors presented their inver- 
sion of the MACHO Alert 95-30 light curves, by seeking a non-negative decreasing solution. 
No stability estimates were made in this work either, moreover the data manipulation is 
questionable (neglecting observational errors, converting fluxes to amplifications prior to 
inversion etc.). 

Even the constrained method has its implicit drawbacks. Either the unconstrained solu- 
tion given by equation (9) automatically satisfies the desired constraints, or the constrained 
solution has to satisfy at least one equality instead of the required inequality. This is a 
direct consequence of the minimized function being quadratic and positive definite in its 
parameters B{rj). Such a function has only one global minimum, given by the unconstrained 
solution (9), and no other local minima. If this minimum does not lie within the constrained 
region of parameter space, the constrained solution lies on the boundary of the region. A 
requirement of positiveness therefore forces the constrained solution to have at least one 
zero point, a requirement of monotonic decrease forces a constant section of the solution, 
concavity forces a linear section, etc. A combination of constraints forces an equality in at 
least one constraint. Even if the method provides reasonable-looking solutions, there is no 
guarantee these approximate the actual brightness profile of the source. Without an error 
analysis based on simulations the results cannot be trusted. 



The Backus- Gilbert method (e.g., Craig & Brown 1986) takes quite a different approach 
to inverting equation (1). It sets out to determine an approximation of the brightness at 
some point Tq as a linear combination of all the observed fiuxes. 



3.2. 



Backus-Gilbert Method 



n 




(11) 



i=l 



^This choice is not ideal, as non-transit events contain little information about the surface brightness 
distribution of the source (see HSL and Appendix A). 
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where //^(ro) are suitable coefficients. Substituting the theoretical expectations F{li) for the 
fluxes from equation (1), we find 

1 

S(ro)~ J S{r,ro)B{r)dr , (12) 



where the introduced spread function S{r, Tq) is a linear combination of the integral kernels 
corresponding to all lens positions in equation (1), using the same coefficients, 

n 

S{r,ro) = '^Hi{ro)A{li,r)r . (13) 
1=1 

The spread function should be normalized so that a uniform brightness profile is recovered 
identically by equation (12), namely S{r,ro)dr — 1. The same equation suggests it is 
desirable to select the coefficients //^(ro) so that the spread function approaches the delta 
function 5{r — ro). This narrowness (high resolution) of the spread function is sought by 
requiring for example 

1 

j{r-rofS^{r,ro)dr = Min. (14) 



The final constraint imposed on the coefficients limits the error magnification, 

n n 

4(.o) = E/^'(^o)^f <M , (15) 

i=l i=l 

the constant M is a measure of the statistical error tolerance. The set of coefficients yU.i(ro) 
can now be easily found analytically by minimizing equation (14) using the definition from 
equation (13) together with the normalization condition and the error magnification con- 
straint. We then compute the brightness point B{ro) from equation (11) and its error due 
to light curve noise from equation (15). 

The method tends to be laborious - the procedure has to be repeated for each desired 
point B{ro) individually. Another drawback lies in having the free parameter M. Decreasing 
its value on one hand increases the stability of the result towards light curve noise, but on 
the other hand decreases the resolution (yields a broader spread function). As a result, 
some user-defined optimization is necessary. It should be noted that there is no explicit 
requirement that the obtained solution provide a good fit to the light curve data. The result 
may be stable towards light curve noise and yet be incorrect due to systematic effects. This 
should be checked independently. 
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Preliminary reports on testing the Backus-Gilbert method were presented by Hendry et 
al. (1998) with a recent update by Gray & Coleman (2000). In both cases their presented 
sample result is for a lens with an Einstein radius e = 1 and a non-transit impact parameter 
p = 1 . It should be noted that for these parameters the surface brightness signature on the 
light curve is very weak (see HSL and Appendix A), in fact lower than the 2% noise added 
to the simulated light curve in Hendry et al. (1998)^°. Using a larger Einstein radius (which 
is also more realistic) should provide better results. At radii 0.6 < r < 0.9 the recovered 
brightness approaches the original profile, though there is a systematic discrepancy in the 
residuals. Between 0.5 and 0.6 the obtained profile peaks and starts decreasing (unphysically) 
towards smaller radii, reaching a ~ —18% deviation at r = 0.05 (even worse, ~ —25% in 
Gray & Coleman 2000). As a consequence, the recovered profile has a lower integrated flux, 
and thus will not provide a good fit to the observed light curve. The authors state that 
the deviation for the region r < 0.6 does not improve with reduced light curve noise. The 
inversion turns unacceptably noisy for 10 or less light curve points according to the authors. 
It is apparent that in this case the actual inversion error is determined by systematic effects 
rather than by statistical errors. The accuracy is too low for the results to be of much use 
for measuring the limb darkening of the source. 

The Backus-Gilbert technique seems to be poorly suited for microlensing light curve 
inversion in general. In the case of non-transit events, all the kernels on the right-hand 
side of equation (13) are very similar functions, increasing monotonically from the source 
center towards the limb. It is difficult to construct a spread function narrowly peaked at any 
selected point on the disk by a linear combination of such terms. At any rate, this would 
require at least some of the coefficients ^{{tq) to be negative. However, subtracting fluxes in 
equation (11) only increases the effects of noise. The problem in recovering the brightness 
profile of the inner disk is also apparent: at r = all the kernels are identically equal to zero. 
Getting a spread function peaked in this region would require high values of the coefficients, 
which would increase its amplitude at the limb even more. 

In the case of transit lens positions, the corresponding kernels in equation (13) are 
divergent at these points. Getting a spread function at least effectively narrow-peaked at a 
particular point by a weighted combination of functions divergent at different points seems 
to be difficult. It remains to be seen whether and how the Backus-Gilbert technique could 
be modified to be useful for microlensing light curve inversion. 



^•^Gray & Coleman (2000) show that the results do not change qualitatively if the light curve noise is 
varied from 0.2% to 10%, which suggests the result is an artifact of the method rather than a measurement 
of the source profile. 
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3.3. Limb Darkening Decomposition 

We can avoid a number of problems encountered by the previous methods by decom- 
posing the brightness profile into a linear combination of nc < n basis functions fj{r), 

B{r)c^B{r) = J2<^jfjir) ■ (16) 
We can now rewrite equation (1), 

F{k) = / B(r) A(k, r)rdr^Yl '^^j'^j ' (1^) 

n .7 = 1 





where the matrix elements Cij are defined by 

1 



Cij= I A{h,r)fj{r)rdr , (18) 
b 

and can be therefore computed to arbitrary precision. We now use equation (17) to minimize 
given by equation (8), obtaining the best- fit coefficients 



n nc 



a, = Y.^{C^C)J,'C,,F, . (19) 



i=l k=l 



As previously, we rescaled the matrix elements to Cij = Qj/cTj. The standard deviations of 
the obtained coefficients due to light curve noise are 

<7«, = ^(C^C)-/ . (20) 



This method is mathematically closely related to inversion by discretization. However, 
unlike in the case of discretization, the involved integrals can be computed to arbitrary pre- 
cision and therefore do not contribute to the inversion error. The quality of this method is 
given by the accuracy of the decomposition in equation (16), which is the only introduced 
approximation. In principle, with a sufficiently high number of decomposition terms the 
range of brightness profile shapes is restricted only negligibly. However, with a limited num- 
ber of terms this method is clearly model-dependent, unlike discretization. In the following 
section we proceed to construct a basis fj{r) suitable for inverting broadband light curves. 
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4. SELECTION OF DECOMPOSITION BASIS 

To test the decomposition method, we generate simulated microlensing hght curves using 
sample brightness profiles obtained directly from model atmosphere computations (see HSL 
and Heyrovsky & Loeb 1997 for more details). We explore a principally four- dimensional 
parameter space, setting the Einstein radius e, impact parameter p (position of the closest 
light curve point here) , number of light curve points n (in most cases only transit points are 
used), and relative noise level 77. We add noise to the simulated light curve by setting 

F, = F(/,)[l + r/,] , (21) 

where 77^ (i = 1, . . . n) is a Gaussian-distributed random number with zero mean and standard 
deviation rj. Hence r]F{li) represents the error cTj. 

For a given basis fj{r),j = l,... ric we evaluate the quality of the inversion {B/B — 1) 
and compare it with results obtained for the standard limb darkening (hereafter SLD) basis 

gi(r)^l, gj(r)^-[l-VT^V^' forj>l. (22) 

In this basis two terms {nc = 2) correspond to the usual linear limb darkening law, three 
to the quadratic law, one to uniform brightness. Four out of the five published microlensing 
limb darkening measurements mentioned in §1 used the linear law. Albrow et al. (1999a) 
measured an additional term (1 — v^l — r^, square-root law) in MACHO 97-BLG-28. 

In the following subsection we illustrate some general features and problems of this 
method on the case of an analytical basis. In §4.2 we derive an ideal basis by principal 
component analysis of a set of brightness profiles obtained from realistic model atmosphere 
calculations. 



4.1. Polynomial Basis 

We first try to set fj{r) = Pj{r), where Pj{r) is a polynomial satisfying the requirement 
dPj/dr{0) ~ 0, -Pi(r) is a constant, and Pj{r) is of degree j for j > l}^ Any brightness 
profile can be well approximated using several terms of such a basis. In the following we 
briefly summarize the main results of inversion tests. 

The brightness profile close to the limb is generally poorly constrained by light curve 
data. For practical purposes we concentrate on the inversion accuracy for the interior region 
r < 0.9 of the stellar disk. 



The inversion result B{r) is independent of the exphcit form of the basis satisfying these criteria. 
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As one would expect from the microlensing chromaticity curve (see Appendix A), the 
lower the impact parameter, the better the inversion. While non-transit points with small 
errors can help improve the inversion accuracy, transit points are necessary for achieving an 
accuracy better than 10%. The inversion quality is independent of the Einstein radius for 
e > 3, it deteriorates slowly for e < 3. This is in agreement with the microlensing sensitivity 
to surface structure, which doesn't increase much for e > 3 while for smaller values^^ it drops 
to zero, as illustrated in Appendix A. 

The inversion accuracy rapidly deteriorates when using more than a couple of basis 
terms. Moreover, the number of terms necessary for obtaining the best inversion is not 
correlated with the changing x^/ (degree of freedom). As one cannot directly judge the 
quality of the inversion by the quality of the light curve fit, it is necessary to predetermine 
an ideal number of terms to extract from a light curve with given parameters. Generally, 
we find it is possible to extract 4 polynomial terms for a noise level of 0.1%, 3 terms for a 
currently achievable level of 0.5%, and merely 2 terms (i.e., a parabola) for more common 
noise levels of 1% and 2%. In the SLD case, only 2 terms (i.e., standard linear limb darkening) 
are reliably obtainable in general. 

Overall, the SLD inversions tend to be superior to the tested polynomial inversions. 

However, the quality of the SLD results depends on the underlying brightness profile B{r). 
The results can be systematically worse in profiles poorly described by linear limb darkening. 
The tested polynomial basis is not ideal cither. It describes a too general class of functions, 
so that it is not possible to extract a sufficient number of terms to provide a reliable and 
useful light curve inversion. From these tests we conclude we need a basis significantly 
more constrained than general polynomials, but capable of describing already in two terms 
a broader range of limb darkening profiles than does standard linear limb darkening. 



4.2. Basis Constructed by Principal Component Analysis 

For constructing an adequate basis we start out with BVRI profiles obtained from 
eight computed red giant model atmospheres (see Appendix B and Fig. 2), using filter 
response functions provided by Bessell (1990). Our goal is to find a basis giving the best 
gradual (single-term, two-term, . . .) approximation to these um — 32 different brightness 
profiles Bi{r), i — 1, . . .tim ■ We can obtain the first normalized component of the basis 



^^As noted in HSL, for practical purposes the lower cutoff is e ~ 0.4. Events with smaller Einstein radii 
are too weak (maximum amplification < 1.34) to trigger an alert in most current surveys. 
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{Jo /i (^) dx — 1) by minimizing the total integrated square residual 

W) = E j Yii^)-[j^iiy)fiy)dy\fi^)] > (23) 



i=i 



where 6j(r) are the model profiles, suitably normalized. The inner integral is the coeffi- 
cient approximating bi{x) by f{x). The minimization can be performed by the calculus of 
variations, leading to an equation of the form 



^ riM 



/"-IVl 
J2k{xMy)f{y)dy = Xff{x) . (24) 



We can see that the sought function /i(r) is an eigenfunction of a positive semidefinite 
integral operator with a symmetric kernel. It can be shown not only that /i(r) corresponds 
to the highest eigenvalue A/^ , but also that all the further components of the basis can be 
obtained by extracting normalized eigenfunctions with successively lower eigenvalues^^. The 
total integrated square residual for a /c-component basis is 

riM \ f k N 2 

Wi,---A) = E J y^i^)-J2[J ^^(y^My)dy]M^)j dx , (25) 

its minimization leads to equation (24) for each of the components independently, thanks to 
their orthogonality. The eigenvalue for the j'-th component is equal to the decrement of the 
residual 31 due to the addition of the j-th component to the basis, 

nM J. 2 

By substituting the constructed eigenfunction basis into equation (25), we obtain the actual 
value of the total residual. 



■■■fk) = Y. ^?(^) d^-Y. ■ (27) 

i=i i j=i 



^^As a bonus, orthogonality of the eigenfunctions, Jg fi{x)fj{x)dx = 6ij , is guaranteed, because the 
operator in equation (24) is Hermitian. 
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Prom this expression we can see that the normahzation of the brightness profiles which gives 
each an equal weight in this analysis is 

1 

h(r)=[jB^{y)dyy'^'B,(r) , (28) 



and consequently , . . . fk) — um — Yl!j=i ^ fj ■ there are nj linearly independent profiles 
among Bi{r), there are exactly ni positive eigenvalues - all the rest are zero. Their sum is 
Um, hence they leave a zero residual. In most cases based on non-trivial model atmosphere 
calculations (such as in this case), all Bi{r) are linearly independent and ni — um- 

The approach sketched above is a special case of principal component analysis (e.g., 
Preisendorfer 1988; JoUiffe 1986). The eigenfunction problem involves solving the homoge- 
neous Predholm integral equation (24), which can be easily achieved numerically by iteration. 
The first three basis functions obtained by principal component analysis of our 32 model pro- 
files are presented in Pigure 3. 

We demonstrate the superiority of the first two components of the basis over linear limb 
darkening in Appendix B, for the model profiles as well as for data from optical observations 
of the Sun. This result suggests that although the basis was constructed from brightness 
profiles of cool red giants, it can be used for a broader class of sources. On the other 
hand, similar bases can be constructed using sets of model atmospheres of different stellar 
types. The main reasons for using red giants in this work are their intrinsic brightness 
(observational advantage of achieving better photometry) and large size (advantage of getting 
better coverage of the transit). In addition, red giants are numerous in the Galactic bulge 
fields observed by the microlensing survey teams. Pinally, the only well-documented single 
lens microlensing transit event by the time of writing, MACHO Alert 95-30 (Alcock et al. 
1997), involved a red giant source. 

Initial inversion tests with the basis obtained by principal component analysis (hereafter 
PCA) and fixed lensing parameters as in the previous subsection show that the systematic 
errors plaguing the SLD inversions are significantly reduced. The obtained results are gen- 
erally more accurate and less biased. However, it is necessary to test PCA inversion in the 
full light curve analysis, without knowing the lensing parameters a priori. 



5. FULL INVERSION OF SIMULATED LIGHT CURVES 



In order to demonstrate the reliability and quality of the described PCA inversion, we 
test the method for a representative range of lensing parameter combinations and source 
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brightness profiles. We estimate the statistical error by inverting a number of noisy light 
curves for each parameter combination, and the systematic error from the results obtained 
for different brightness profiles. For comparison, we perform the same procedure using linear 
limb darkening (hereafter LLD). 

In the following subsection we describe the implementation of the simulations and the 
light curve inversion algorithm. In §5.2 wc discuss the accuracy of the recovered lensing 
parameters, and in §5.3 the obtained limb darkening of the source. Based on these results, 
we suggest in §5.4 how to analyze observed light curves in order to accurately recover the 
limb darkening as well as the lensing parameters. 



5.1. Simulation and Inversion Algorithms 

We generate simulated light curves for combinations of five lensing parameters: Einstein 
radius e, impact parameter p, blend parameter f3 (ratio of flux from unresolved stars to 
baseline fiux F*(oo) of lensed star), number of points per source diameter crossing time 
(equal to transit point number for small impact parameters), and light curve noise level 
f]. The two parameters defining the time scale of the simulated event can be kept fixed at 
arbitrary values, for simplicity we set the source radius crossing time U — 1 and the time 
of closest approach to — 0. However, we do fit for both of these parameters in the following 
light curve analysis. In order to get a reasonable fit for most of the parameters we include 
non-transit points, setting the total number of points to n = 2n^, and distributing them 
in equal time intervals from roughly 2t^ before the peak to 2t^ after the peak^^. The lens 
positions are given by 

k = Vp' + [{ti-to)/U]^ (29) 

for i = 1, . . . n. We compute the source star fluxes F^{li) directly from equation (1) using one 
of the model brightness profiles (not its approximation) and the angle-integrated amplifica- 
tion from equation (4). We integrate across the transit point divergence using expression (6) 
and, to match the accuracy of this expansion, a local linear approximation of the brightness 
profile. To obtain the simulated fiuxes Fj we add the blend contribution and Gaussian noise 
Tji with standard deviation rj (as in the previous section), 

Fi^[F,{k) + pF,{oo)][l + 'ni] . (30) 

The errors are then given by (jj — r] [F^{li) + P F^{oo)]. 



^^Explicitly, ti = to + [(4i — 3)/n — 2] , i = 1, . . . n . This distribution is slightly asymmetric around the 
peak. 
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The studied parameter space is a Cartesian grid defined by all combinations of values 
e E {1,2,5,10}; p E {0,0.25,0.5,0.75}; f3 E {0,0.25,0.5}; n, E {10,20,40}; and E 
{0.005, 0.01, 0.02}. As shown in Appendix A, higher values of e provide the same sensitivity 
to surface structure as e = 10, even though they provide higher overall amplification. For 
lenses with Einstein radii lower than 1 the sensitivity drops, for even lower values such 
events would be hard to detect at all, due to their low amplification. Events with impact 
parameters p > 1 contain only a weak signature of the surface structure, so that their 
inversion is useless under realistic noise conditions, as shown in the previous section. The 
choice of blend parameters is geared toward events with bright (e.g., red giant) source stars, 
for which the fraction of blended fight can be expected to be low. The point numbers are 
representative and realistic. The 0.5% noise level has been achieved in several microlensing 
events, the two other values are more usual. 

The model fitted to the fight curves is based on equation (17) with additional blended 
fiux Fb, 

Fuih) ^ FB + OiiCii{e, li) + a2Ci2{e, k) , (31) 

where the matrix elements Cij are given by equation (18), and the lens positions by equa- 
tion (29). In this model the asymptotic flux (flux in absence of the lens) from the lensed 

star is -F*(oo) = ai Coi + ol2 C02, where Coj = lim^^^oo Cij for j = 1, 2 . These two values are 
Coi = 2.680, C02 = 2.418 for the PCA basis and Cqi = 3.142, C02 = -1.047 for the LLD 
basis. The value Fb/F*(oo) corresponds to the blend parameter /3 as defined above. 

This seven-parameter model is linear in three parameters {Fb, «!, 02), and depends 
on the remaining four (e, p, t^, to) through Cjj[e, k {p, t*, to)]- In order to speed up the 
subsequent computations, it is crucial to precalculate these matrix elements for a sufficiently 
dense grid of Einstein radii e and lens positions / using the two basis functions (/i 2 for PCA 
and gi^2 for LLD) and equations (4) and (6), in a manner similar to the source light curve 
computation described above. Values at intermediate points can be obtained by interpola- 
tion. 

Due to the form of the dependence on the nonlinear parameters the full inversion prob- 
lem cannot be solved analytically. Instead we find numerically the best-fit parameters by 
minimizing 

n 

x' = J2<'^"i^^-^M{k)r ■ (32) 

i=l 

However, the numerical search can be limited to only four out of the seven dimensions, due 
to the presence of three finear parameters (as also noted by Rhie et al. 1999 and Albrow 
et al. 1999b). For any fixed combination of the four nonlinear parameters, the combination 
of (Fb, CKi, a2) minimizing the function can be found analytically by solving a set of 
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three linear equations. These are the corresponding minimum conditions obtained by dif- 
ferentiating equation (32). It is advantageous to exclude non-physical solutions in advance 
by requiring for example a non-negative blend, a non-negative brightness profile peaked at 
the center (see Appendix B), and the total model asymptotic flux being within ±o"o of the 
baseline flux Fq (both constants estimated from the data). These flve constraints restrict 
the linear parameter space to a pentahedron. If the computed minimum (for any fixed set 
of nonlinear parameters) lies outside this region, the constrained minimum can be found on 
its surface. 

For the minimum search in the remaining four-dimensional parameter space (e, p, t*, to) 
we used a Monte Carlo technique based on gradual contraction of parameter intervals towards 
lower values of x^- This technique was found to be very robust in finding the global minimum. 

As the full inversion problem is nonlinear, we inverted 100 noisy light curves for each 
point of the parameter grid to get a realistic estimate of statistical errors. In order to find 
the systematic error and bias of the method we tested the 32 model brightness profiles fisted 
in Appendix B. The distribution of the values recovered from different noisy simulations 
was usually found to be substantially non-Gaussian. Therefore we decided to describe these 
distributions using the median instead of the mean, and the narrowest interval containing 
68.3% of the recovered values instead of the ±0" standard deviation interval, which contains 
the same fraction of points in a Gaussian distribution. 

With respect to the two photometry-defining parameters rj and n* , a good-quality fit 
not dominated by systematic errors is indicated also by the linear scaling of the total error 

— 1/2 

with light curve noise and the inverse square root of the number of points (oc T^n* ). We 
tested both these scalings; for the scaling we used the results for rj — 0.005 in order to 
suppress the effect of light curve noise. 

5.2. Recovered Lensing Parameters 

Before studying the recovered limb darkening, we check here the accuracy of recov- 
ered lensing parameters: namely, Einstein radius e, impact parameter p, blend parameter 
/3 = -Fb/(«iCoi + 0:2^02)1 source radius crossing time and closest approach time to- 
For avoiding confusion we add a subscript '0' to denote the original values used for the 
simulations, i.e., co , Po > /^o > and too- For each point in the original parameter grid, we 
statistically analyze the full set of 32 x 100 recovered parameter values. 



These two descriptions converge in case the underlying distribution is indeed Gaussian. 
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We find that the best results for most parameters are obtained in events with high 
Einstein radius, low impact parameter and low blending. For example, the original parameter 
combination {eo = 10; po = 0; /Sq = 0;t*o = l;^oo = 0} was recovered using the PCA basis 
as {e = 10.11+°:}^;^ = 0.016+^16;/?= (l-8l?i) x 10-';U = 0.9989^^^61; = (-O-ltH) x 
10~^} from simulated light curves with n* = 40 and rj = 0.01. As an exception, the impact 
parameter is recovered most accurately in events with higher impact parameters, transiting 
closer to the limb of the source star. This is probably due to the weak dependence of 
the microlensing amphfication on the brightness profile in the range 0.65 ^ I O.S (see 
Appendices A and B), which in turn implies a lower degeneracy between the remaining 
fitted parameters (mainly p and ) . 

The closest approach time to is by far the best-recovered parameter. The accuracy of 
the recovered Einstein radius and blend parameter decreases slightly in events with eo ^ 5. 
In such events, better accuracy is obtained by including more distant light curve points in 

— 1/2 

the analysis, as discussed further in §5.4. The total error of the PCA results is oc jyn* 
in events with eo > 1. The only found exceptions are the impact parameter error when 
Po = 0, and the blend parameter error when Pq — (both are values at the parameter space 
boundary) . 

The LLD results are in most cases qualitatively similar. In a few cases the accuracy is 
slightly worse and the error scaling with light curve noise and number of points docs not 

— 1/2 

follow r^n* as accurately as in the PCA analysis. Clear systematic diflcrcnccs can be 
seen in the recovery of the source radius crossing time t^,. The most interesting difference 
is found in the general scaling of the errorbars: the errors scale roughly as 7]^'^ and n^°'^, 
quite consistently. A closer study of the numerical results reveals that systematic errors in 
the values recovered using LLD are large - comparable to the statistical errors, hence 
the different scalings. This means that while the LLD errors for — 40 and rj — 0.01 are 
only slightly larger than the PCA errors, decreasing the light curve noise or increasing the 
number of points will only boost the relative advantage of PCA analysis over LLD analysis. 
This problem is most probably caused by LLD generally providing a poor fit to realistic 
brightness profiles, as demonstrated in Appendix B. 



5.3. Recovered Limb Darkening 

In this subsection we first present the results obtained for three sample input brightness 
profiles. We then proceed to estimate the actual error of the reconstructed profile for a 
source with unknown limb darkening by a combined analysis of the simulation results for all 
32 tested profiles. 
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The shape of the recovered profile is given by the PCA or LLD hmb darkening parameter 
K = a2 /cKi . The accuracy of recovering this parameter can be studied similarly as the 
lensing parameters in the previous subsection. Of the 32 input model brightness profiles (see 
Appendix B and Table 1), we present in Figure 4 the results for the most centrally peaked 
profile (B3750gl0), an intermediate one (R3500g05), and the flattest proflle (MOOOglO). We 
determine the median and the 68.3% confidence interval from the 100 noisy light curves 
computed for each input parameter combination. In addition we plot the minimum and 
maximum obtained values (none of the simulation results were discarded). The results 
shown in the figure correspond to n* = 40 and rj — 0.01 . 

As seen from the results for all profiles and both PCA and LLD fits, limb darkening 
parameters are most accurately obtained for events with low impact parameter pq, although 
the dependence is weak for po < 0.5. The accuracy increases only very weakly with higher 
Einstein radius eo, and is fully independent of blend parameter /3o. These results are en- 
couraging - either limb darkening parameter can be recovered even if the blend or Einstein 
radius are poorly fitted. 

In the case of profile B3750gl0, the original LLD parameter value (1.0) hes exactly 
at the parameter space boundary, as refiected also by the recovered values. The plotted 

PCA errorbars for the ideal lensing parameter combination (eo = 10 ; po = ; /3o = 0) are 
consistent with best-fit PCA limb darkening parameters corresponding to all the B profiles 
in Table 1. In the LLD case, the errorbars arc formally consistent only with the B3750gl0 
value, due to their nearly zero width. However, this does not mean that the actual profile is 
recovered better (sec Table 1). 

The results for profile R3500g05 exhibit a nonlinear dependence of the LLD errors on the 
light curve noise rj, indicating non- negligible systematic errors. The LLD medians are clearly 
correlated with the impact parameter Pq. The PCA errorbars for the ideal lensing parameter 
combination are consistent with best-fit PCA limb darkening parameters corresponding to 
the R3500 and V4000 profiles in Table 1. In the LLD case, the errorbars are consistent with 
all profiles from R3500 to B3500, including all V profiles. 

The recovered LLD medians for profile MOOOglO systematically underestimate the orig- 
inal value. The PCA errorbars for the ideal lensing parameter combination are consistent 
with best-fit PCA limb darkening parameters corresponding only to the 14000 profiles in 
Table 1. The same result is obtained in the LLD case. However, note here that no flatter 
proflles were used (see Fig. 2 and Table 1). 

From the general analysis we conclude that the PCA results have the advantage of being 
uncorrelated with the lensing parameters and having no apparent systematic problems. The 
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LLD medians are correlated with po and sometimes with /3o ; the results suffer from systematic 
errors due to the inadequacy of the LLD approximation. Even in ideal cases, the LLD results 
are less accurate than the PCA results. 

It is worth noting that this test only told us how well a fitted parameter of the actual 
brightness profile can be recovered from the light curve. The conclusions confirm that the 
parameter of a worse-fitting model will be recovered worse. The following part of this 
subsection studies how well is the actual brightness profile of the source recovered. 

Using the results for all 32 profiles it is now possible to make an estimate of the accuracy 
of the inverted brightness profile B(r), given by equation (16), for any given combination 
of input parameters (eo , po , f3o , , r]) independent of the actual source brightness profile. 
First we compute the relative brightness profile deviation for each of the 32x100 inverted 
light curves, 

F,o(oo) B{r) i^*o(oo) ai/i(r) + a2/2(r) ' ^ ^ 

here F*o(oo) is the asymptotic flux of the input model profile B{r) Prom this full set of 
data we determine the total median 5s(r), which represents the bias of the inversion, and 
the 68.3% confidence interval {Sb — Sl, 5b + Su) , where + Su{f) is the total error. 

To estimate the systematic error, we determine the median deviation for each model 
profile separately, from its 100 noisy light curves. Analysis of the distribution of these 32 
profile-specific medians yields a "median median" 5mm('"), and the 68.3% confidence interval 
i^MM — ^SL , ^MM + ^su) ■ The value Ssl{t^) + SsuIt^) is the systematic error, a measure of 
source-dependence of the inversion accuracy. 

Prom the previous analysis we obtain a best estimate of the brightness profile of the 
lensed star, 

+Suir) 



B(r) 



l + 5B(r) 



X B{r) . (34) 



If the accuracy of the recovered blend parameter is low, by normalizing this result we still 
obtain the best estimate of the normalized profile of the star. 

A reliable inversion method requires not only a small total error Sl + Su, but also 
neghgible bias Sb/{Sl + Su) — > 0, so that the inversion result is closest to the best estimate. 
It is also desirable to have small systematic errors {Ssl + Ssu)/{Sl + Su) <S 1 , so that the 



^^The obtained profiles B{r) are normalized to the total source flux in order to avoid uncertainties intro- 
duced by any potential inaccuracy of the fitted blend. 
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inversion accuracy doesn't depend substantially on the source profile, and the inversion can 
yield better results with improved photometry. Such a method would also have 5mm — Sb- 

Representative PCA and LLD results of the described analysis are presented in Figure 5 
for eo = 10, Po — 0, n^ — 40 and rj — 0.01. The plotted curves include the bias Sb, the total 
error as the 68.3% confidence interval {Sb — Sl, Sb + Su) , the "median median" 5mm and the 
systematic error as the 68.3% confidence interval {Smm — Ssl , 5mm + 5su) ; all are plotted 
as a function of the position on the stellar disk r. The jagged shape of some of the curves 
is merely due to using only 32 model profiles in the analysis. All the desirable properties 
mentioned in the previous paragraph can be achieved by PCA inversion, but not by LLD 
inversion. The PCA results show a negligible bias all the way to the limb of the source, 
the coincidence of medians 5mm — 5b, and systematic errors much smaller than total errors 
- with a single exception in the case of the best-fit point at r ~ 0.7 with an error on the 
order of 0.5% The quality of the inversion improves weakly with higher Einstein radius 
eo, and is practically independent of the blend (Sq (only in poorly fitted events the stability 
deteriorates very weakly with increasing /5o). The errors scale linearly with light curve noise 
T) and decrease as the inverse square root of the number of points n^, per source diameter 
crossing time. The best inversion is achieved for low impact parameter events. However, 
even at po = 0.75 good accuracy can be achieved by obtaining a higher number of good 
quahty transit points. The results in Figure 5 demonstrate the best obtained inversion (for 
Po = 0) has errors < 2% from the disk center out to r ~ 0.9. These errors can be further 
reduced by having better than 1% photometry or more than 40 light curve points per source 
diameter crossing time. 

While the LLD inversion results are not bad either, they do not achieve the quality of 
the PCA results and do not have the desirable inversion properties. The results are always 
biased, on average producing flatter profiles. The bias can reach as much as 2% in the region 
r < 0.9 and diverges further closer to the limb. The systematic errors are large and often 
dominate the total error, particularly in the region 0.5 ^ r < 0.8 . The two medians 5mm 
and 5b differ by as much as 1%. The trends of the results with cq, po and f3o are the same as 
in the PCA case. The linear scaling of the errors with rj and their inverse square root scaling 
with occurs in fewer cases and to a lesser extent than in PCA inversion. The quality of 
the achievable results is limited by the relatively large systematic errors. 



^''This point occurs in the results merely due to the fact that the intersection point of the PCA limb 
darkening model (as discussed in Appendices B and A) lies in the narrow intersection region of the used set 
of model profiles. This is not the case for the LLD model. 
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5.4. Inversion of Observed Light Curves 

In the simulated inversions performed in the previous subsections we only used points 
near the light curve peak, within ±2 source radius crossing times from the closest approach. 
Light curves of actual niicrolensing events usually provide data going well beyond this region. 
To study the effect of more distant points on the recovery of lensing parameters and the limb 
darkening, we performed simulations including light curve points within ±2 Einstein radius 
crossing times from the closest approach. 

The results for eo ^ 5 show a dramatic improvement in the recovered Einstein radii and 
blend parameters, as mentioned above. In such events impact parameters are recovered bet- 
ter as well, with the exception of the already well recovered po = 0.75 which improves only 
marginally. The distributions of the recovered source radius crossing times have more accu- 
rate medians, but the total errorbars remain unchanged. On the other hand, the accuracy 
of the recovered closest approach time and (more notably) the limb darkening parameter is 
worse than while fitting only the event peak. Consequently, the total error of the recovered 
brightness profile at source center increases by 1-3%. 

These findings are a natural result of the light curve fitting procedure. While fitting an 
extended light curve going far beyond the transit region, the of the fit is altered more by 
a minute change in the global parameters (e.g., e and j3) than by a more significant change 
in the limb darkening parameter, which affects only points in the transit region - a small 
fraction of the total number of light curve points. 

High accuracy for the limb darkening of the source as well as the lensing parameters 
can be achieved by a two-step procedure. In the first step we fit the full light curve, varying 
all parameters. In the second we fit only the transit region (not more than ±2 source 
radius crossing times from the peak), varying the limb darkening parameter and the closest 
approach time, keeping the other parameters fixed at values obtained in the previous step. 
Finally, we check the value of the overall result. 

6. ANALYSIS OF MACHO ALERT 95-30 LIGHT CURVES 

In this section we analyze the observed light curves of the microlensing source-transit 
event MACHO Alert 95-30 (Alcock et al. 1997). Altogether seven light curves were measured 
by the MACHO and Global Microlensing Alert Network (OMAN; Becker et al. 1997) collab- 
orations. The individual hght curves were obtained by the MACHO team at Mount Stromlo 
Observatory (MSO; MACHO B filter - 293 points used, MACHO R filter - 205 points), and 
by the GMAN collaboration at Cerro Tololo Inter-American Observatory (hereafter CTIO; 



-24- 



R - 106 points), Mount John University Observatory (MJUO; R - 41 points), University of 
Toronto Southern Observatory (UTSO; R- 133 points, V- 127 points) and Wise Observatory 
(WISE; R - 17 points). Alcock et al. (1997) presented the results of hght curve fits using a 
uniform source model as well as a linearly limb-darkened source with fixed limb darkening 
parameters based on theoretical expectations. Here we redo the general light curve fitting 
procedure^^ to measure the limb darkening profile of the source, as well as to obtain more 
accurate lensing parameters. 

The analysis can be performed in the manner used in the previous section. The only 
difference is that this case involves a simultaneous inversion of seven light curves. The four 
nonlinear lensing parameters (e, p, t^, to) are common to all the curves, the linear param- 
eters (Fb, cti, 0:2) should be determined for each light curve separately. This is necessary 
because of differences not only in the used filters, but also in the optics, CCDs, point spread 
functions, data reduction methods etc. Fitting the 25 parameters to the full set of data is 
straightforward. As described in the previous section, for a given combination of the four 
nonlinear parameters all the linear parameters can be computed analytically by minimizing 
the total function. Errors of the obtained best-fit parameters can be obtained by Monte 
Carlo simulations. We generate 100 sets of seven noisy light curves using the obtained best- 
fit parameters, plus the timings and errorbars of the original data. We invert these curves 
and from the distribution of the obtained values we determine the (non-Gaussian) 68.3% 
confidence intervals for each parameter. 

Estimation of the baseline parameters Fq and cxq for the two MACHO light curves 
revealed a systematic problem. As shown in Figure 6, the baselines exhibit nearly periodic 
correlated variations. There is an apparent period of ~ 45 days, with an amplitude of 
~ 3.5% in MACHO B and ~ 2.9% in MACHO R. Further inspection of the hght curves 
suggests there is additional irregular variability on both longer and shorter timescales. Due 
to this complex variability, even within the span of 10 days around the light curve peak the 
intrinsic fiux from the source may change by ~ 3%. As an immediate consequence, this 
change will obscure the weak chromatic lensing effect close to the source limb (see Appendix 
A). Moreover, the limb darkening parameters (i.e., the source brightness profiles) can be 
expected to vary, rendering any inversion questionable in principle. Finally, the obtained 
lensing parameters will also be inaccurate to a certain degree, due to the source variations 
in the light curve peak region. 

Proceeding formally with the fit, the problems become even more apparent. The 



^^The event data were kindly provided by Andy Becker, who re-reduced the CTIO, UTSO and WISE hght 
curves recently. 
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value for the fit to the original data is substantially and systematically higher than the 
values obtained for Monte Carlo simulated light curves with the best-fit parameters. Even 
the residuals of the fit close to the peak show systematic deviations. The relative residuals 
Fi/FM{h) ~ 1 are plotted in Figure 7 for this part of the light curve. As can be seen, the fairly 
large residuals in the different light curves are clearly correlated (e.g., note the post-peak 
changes) . 

The formal best fit parameters are: Einstein radius e = 13.9llg o| ; impact parameter 
p = 0.7253^QQQ44 ; source radius crossing time = 2.5164]'^QgQ28 ^ closest approach 

time to = 1321.2359lo;oo38 d, measured from JD 2,448,623.50 (Alcock ct al. 1997). Except 
for to, these values are slightly different from those obtained by Alcock et al. (1997). Partly 
this is caused by the greater number of points near the light curve peak in the current set of 
data (mainly in the UTSO hght curves) , partly by the more general fitting procedure. The 
greatest effect, however, is due to source variabihty. Neither of the fitted models accounts for 
the effect of an intrinsically irregularly variable source. All the lensing parameters obtained 
here as well as in Alcock et al. (1997) should therefore be taken with a grain of salt. 

The values of the formal best-fit limb darkening parameters are concentrated at the 
flat-profile boundary of the explored parameter space. This can be seen also from the largely 
positive residuals during the transit. Moreover, the best-fit limb darkening parameters are 
mutually physically inconsistent (implying flatter profllcs in V than in R etc.). All these 
features conflrm the conclusion that limb darkening measurement from the MACHO Alert 
95-30 event is not possible due to the variability of the Icnscd star. It cannot be ruled out 
that some of the residuals close to the peak may be due to the lens being binary (Griest 
& Safizadeh 1998; Gaudi & Sackett 2000), source being binary (Griest & Hu 1992; Han & 
Gould 1997; Han & Jeong 1998) or even due to asymmetric surface features on the source 
star (Heyrovsky & Sasselov 2000; Bryce & Hendry 2000). However, attempts to fit any of 
these models will also be largely hindered by the variabihty of the source star. 

It should come as no surprise that the lensed M4 giant source in MACHO Alert 95-30 is 
variable. Percy & Polano (1998) state that "virtually every M giant is variable in brightness" , 
with small-amplitude red variables being two orders of magnitude more common than Miras. 
Variability larger than 1% in the V band sets in already in late K and early M type giants, and 
increases further toward late M type giants (Jorissen et al. 1997; Percy & Parkcs 1998). The 
variability seems to be associated primarily with pulsations (Percy & Parkes 1998; Jorissen 
et al. 1997), as suggested also by small-amplitude radial velocity variations measured by 
Cummings et al. (1998). The variations are not strictly periodic, sometimes they are multi- 
modal or irregular (Percy & Polano 1998). The periods can be on the order of tens of days 
down to single days, as recently demonstrated from Hipparcos data by Koen & Laney (2000). 
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These results demonstrate the importance of checking the basehne of source-transit 
event hght curves to detect even small-amphtude variabihty. Only events in which the 
change due to intrinsic variability over the event timescale is sufficiently smaller than the 
expected microlensing chromaticity (see Appendix A) can be used for measuring the limb 
darkening of the source star. As an exception, if the pulsations are well defined, regular and 
stable, they can be accounted for a priori in the lensing model, and the lensing effect can be 
fully separated from the intrinsic variability. MACHO Alert 95-30 unfortunately does not 
fall into this category of events. 

7. DISCUSSION 

The accuracy of measuring the limb darkening parameter from light curves of eclipsing 
binaries was studied in detail by Popper (1984). On the basis of light curve simulations 
Popper concluded that achieving an absolute accuracy of ±0.1 in the linear limb darkening 
parameter requires on the order of 100 light curve points within light curve minima obtained 
with 0.5% photometry. As shown in Figure 4, analysis of microlensing transit light curves 
frequently achieves better accuracy even when using LLD inversion of light curves measured 
with 1% photometry. These results are surpassed in quality by PCA inversion, which in fact 
directly demonstrates the inadequacy of the linear limb darkening model. Going one step 
further, the results in Figure 5 show that PCA inversion can retrieve realistic brightness 
profiles (instead of profile fits) with an accuracy higher than 2% from the center of the 
stellar disk out to r = 0.9 . Inversion of microlensing light curves thus gives the opportunity 
to proceed from measuring the limb darkening parameter to measuring the actual limb 
darkening profile of the source, thanks to the high achievable accuracy and low systematic 
errors. 

Eclipsing binaries have the advantage of having periodic light curves, giving thus fre- 
quent opportunities to measure the eclipse minima. Even though microlensing transit events 
are non-repeating, giving only a single chance to obtain high quality photometry, their anal- 
ysis can place much better constraints on the limb darkening as well as other properties 
of the source (HSL, Heyrovsky & Sasselov 2000). Such events can also potentially provide 
results for a broad range of stellar sources in the currently surveyed Galactic bulge or Mag- 
ellanic Cloud fields, as microlensing does not depend on the stellar type of the source. The 
only biases are introduced by the requirement of good photometry (giving an advantage to 
brighter sources) and a sufficient number of transit points (giving an advantage to longer 
transits, i.e., larger sources). 

The PCA method described in this work on the example of broadband profiles of red 
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giants can be simply extended to any other class of sources, narrow-band filters etc. Using a 
set of brightness profiles obtained from a corresponding class of atmosphere models for the 
required wavelengths or passbands, an adequate PCA basis can be constructed as described 
in §4.2. The first few terms describing a sufficiently broad range of possible brightness profiles 
can then be used for the light curve analysis. 

Observed events with high quality photometry can also be used to discriminate between 
classes of stellar atmosphere models with different PCA bases. The more suitable class 
can be distinguished mainly by the pattern of residuals of the light curve fit. The best-fit 
model should then be checked whether it is in agreement with measiued properties of the 
unamplified source - such as colors, fine strengths etc. Microlensing thus provides unique 
new tests that will contribute to our understanding of stellar atmospheres. 

Microlensing transit events are much more frequent for caustics of binary lenses than 
for single point lenses studied in this work. The reason is that the probability of a caustic 
transit of the source is usually given by the angular size of the caustic rather than by the 
angular size of the source. While the source size still determines the duration of the transit, 
detection of the source entering a caustic gives an opportunity to predict and observe in 
detail at least the exit from the caustic. Limb darkening parameters have already been 
measured for four caustic-crossing events in this way, three of them by The Probing Lensing 
Anomalies NETwork (PLANET) collaboration and one by a joint effort of five microlensing 
collaborations (Afonso et al. 2000). The most recent one at the time of writing was described 
by Albrow et al. (2001). The sensitivity of a linear caustic-crossing event to the brightness 
profile of the source was recently studied by Gaudi & Gould (1999) and Rhie & Bennett 
(1999). As shown by Rhie & Bennett (1999), the chromaticity is generally lower than for 
a single lens. On the other hand, a linear caustic always transits the source completely 
(with "zero impact parameter"). While single lens transits have the potential to deliver 
even better results, binary caustic-crossing events are already delivering good results. The 
approach demonstrated here for the single lens case can be readily applied to the binary case. 
The linearity of the PCA approach will be particularly beneficial, as the fitting of binary 
microlensing light curves is computationally heavily burdened by proliferating nonlinear 
parameters (Albrow et al. 1999b, 2000). This method will substantially improve the analysis 
of all previously reported as well as upcoming caustic-crossing microlensing events. 

8. SUMMARY 

In this paper we show that limb darkening profiles of source stars in microlensing transit 
events can be measured by light curve inversion using brightness profile decomposition. Other 
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inversion techniques at present fail to reach the accuracy achieved by this method. 

Initial tests in §4.1 show that under the usual photometric accuracy only two terms of 
the decomposition can be reasonably recovered. Using linear limb darkening, which suffers 
from systematic errors, is not ideal. It is necessary to find a limb darkening basis providing in 
two terms a better description of realistic brightness profiles than does linear limb darkening 
(LLD). Such a basis can be constructed by principal component analysis of a set of realistic 
model brightness profiles, as demonstrated in §4.2 . The obtained PCA basis provides the 
best possible linear description of the used set of profiles, as demonstrated in Appendix B. 
This method of describing limb darkening should in fact replace the ill-suited and obsolete 
standard power-law limb darkening model in other practical applications as well. 

We demonstrate the described inversion method in §5 by fully analyzing simulated light 
curves. In §5.2 we check the accuracy of recovering the basic lensing parameters. We show 
that best results are obtained for events with Einstein radii e > 5, low impact parameters 
p < 0.5 and negligible blend parameters /5 -C 1. As demonstrated in §5.3, the limb darkening 
parameter is recovered more reliably using the PCA basis than using the LLD basis, achieving 
also better brightness profile resolution. The quality of these results is further demonstrated 
by directly measuring the inversion accuracy, the difference between the recovered and the 
original brightness profile. The best results are obtained for small impact parameters p < 0.5. 
They are only very weakly dependent on the Einstein radius, and independent of the blend. 
For events with small impact parameters even light curves with 1% photometry can yield 
an inversion accuracy of ±2% from the stellar disk center to r = 0.9 near the limb using 
the PCA basis. The accuracy can be further improved by better photometry as well as by 
increasing the number of measured transit light curve points. 

The measured light curves of the microlensing transit event MACHO Alert 95-30 (Alcock 
et al. 1997) reveal the intrinsic variability of the lensed red giant, as demonstrated in §6. 
The time scale, amplitude and irregularity of the variability indicate that the source changes 
significantly at the light curve peak, and thus its limb darkening cannot be well constrained. 
This example underscores the need to check the baseline of microlensing transit events for 
signs of variability as low as 1% on the event timescale. Inversion results cannot be considered 
rehable without this safety test. 

Binary microlensing events, in which a caustic transits the source star, are much more 
frequent than the single lens events studied here. Though their sensitivity to the brightness 
profile of the source tends to be lower, such events have already been used to measure the 
limb darkening parameters of several stars. The analysis of these caustic-crossing events will 
be significantly improved by using the PCA approach to limb darkening as described in this 
work. 
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A. BROADBAND LIGHT CURVE CHROMATICITY 

Following HSL, Figure 8 demonstrates the dependence of the microlensing amplification 
on the surface brightness profile of the source. In this example a red giant source (effective 
temperature T=3750 K, surface gravity logg = 0.5) is transited by an e = 10 lens at zero 
impact parameter p. BVRI light curves of this event in the lower panel are to be compared 
with the corresponding normalized brightness profiles in the upper panel (see Appendix B 
for details). The microlensing amphfication clearly traces the shape of the brightness profile. 
More centrally peaked profiles have more centrally peaked amplification curves (light curves 
for p — 0). Profiles which have a relatively higher amplification at the source center mostly 
have a relatively lower amplification close to the source limb, and vice versa. 

The sensitivity of microlensing to the brightness profile of the source can be measured by 
the relative variation of amplification with wavelength, due to the wavelength dependence of 
the profiles. Following HSL, we define the microlensing chromaticity as the ratio of the total 
amplification range to the average amplification [Ama^. — A„iin]/^aver at a given lens position. 
The maximum, minimum and average amplifications arc evaluated for a particular source 
over a given spectral range at a given spectral resolution. HSL presented the chromaticity 
of low and high resolution spectra, here we demonstrate the broadband chromaticity using 
merely the four profiles from Figure 8. Plotted as a function of lens position in Figure 9, 
the chromaticity curve has the same character as those presented in HSL with a slightly 
lower amplitude due to the lower resolution. The individual curves correspond to lenses 
with different Einstein radii. Note that the chromaticity becomes independent of e for 
e > 5 . The highest sensitivity to limb darkening occurs when the lens is close to the source 
center, another somewhat less sensitive region is close to the source limb. The sensitivity 
drops rapidly for non-transit lens positions. These overall features were first reported in the 
works of Witt (1995); Simmons, Newsam, & Willis (1995); Valls-Gabaud (1995); Bogdanov 
& Cherepashchuk (1995) and Gould & Welch (1996). 

In the vicinity of / ~ 0.73 the amplification is largely achromatic, i.e., independent of 
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the underlying limb darkening profile. The occurrence of this achromatic point is a direct 
consequence of the narrowly confined region around r ~ 0.69 where all the normalized 
brightness profiles in the upper panel of Figure 8 intersect. As Figure 2 and even the higher 
resolution results in HSL demonstrate, this narrow intersection region is common to all the 
computed profiles (see also Appendix B). It can be shown that this property is an indication 
that the full range of profiles is well modelled by a linear combination of two radial functions. 
Any such a suitable model has an exact intersection point common to all the profiles. The 
position of this point, however, depends on the particular model. As a consequence, using a 
limb darkening model with a wrong intersection point in the analysis of microlensing light 
curves will additionally result in inaccurate values of the inferred microlensing parameters. 
However, as mentioned in §7 such a wrong model could be told by the pattern of residuals 
of the light-curve fit. See Appendix B for more details on limb darkening profiles. 



B. LIMB DARKENING PROFILES AND THEIR FITS 

In this work we used limb darkening profiles obtained from eight model atmospheres of 
cool red giants computed by Dimitar Sasselov (see HSL for details) . The effective tempera- 
tures and surface gravities of these models are Te {3500 K, 3750 K, 4000 K} and logg e {0, 
0.5, 1.0}. Specific intensity was computed for 17 rays passing through the atmosphere at 
different positions on the projected disk. These positions arc marked in Figure 2. The 
BVRI profiles were obtained from these results using the filter response functions provided 
by Bessell (1990). To obtain a full radial profile, these 17 points were interpolated using a 
cubic spline with boundary conditions -B'(O) = and B"{1) — The variation of the 
profiles in the studied physical range is strongest with wavelength (filter) , weaker with tem- 
perature and weakest with surface gravity. All profiles with log ^^=0.5 are plotted in Figure 2, 
normalized to unit total fiux. Two trends are obvious: the profiles at longer wavelengths are 
fiatter; so are mostly the profiles of higher temperature models in the studied range. 

The two hmb darkening models used in this paper (PCA, LLD) are demonstrated in 
Figure 10. The PCA basis was obtained directly by principal component analysis of all 32 
computed model profiles, as described in §4.2. The standard linear limb darkening (LLD) 
basis is defined by equation (22). Figure 10 compares the full ranges of profiles covered by the 
two models, normalized to unit flux and plotted on the same scale as Figure 2. The ranges 
are limited physically by requiring positiveness and a central peak for the profile. Notice 



^^The latter condition is fairly arbitrary, but has no global effect on the quality of the interpolation (see 
Fig. 2) due to the high number of computed points constraining the profile close to the limb. 
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that the LLD model is particularly ill-suited for describing the peaked B profiles in Figure 2, 
as well as the behavior close to the limb in general. Also note that LLD profiles with low 
values of the parameter k = a2 /ai are unrealistically fiat. Finally, there is a difference in 
the position of the intersection point as discussed in Appendix A. The LLD point is > 0.05 
closer to the limb than the PCA point, which lies directly in the narrow intersection region 
of the profiles in Figure 2. As a result, simulated microlensing light curves generated using 
LLD will intersect at I ~ 0.77, those generated using PCA at I ~ 0.73 . This is the probable 
cause for the large systematic error in the source radius crossing time recovered from 
microlensing light curves using LLD in §5.2. 

The difference in the quality of least-squares fits^° to all the 32 computed profiles using 
the two limb darkening models is demonstrated by the results in Table 1. The columns 
of the table contain for both limb darkening models the best-fit limb darkening parameter 
K — a2 /oii , the normalized root-mean-square deviation of the fit, 

5fu = {/[S(r) - ai/i(r) - a2f2ir)Ydr / f B\r)drY/' , 
and the relative flux excess AF/F of the best-fit solution over the original profile. The PCA 
fits always leave an r.m.s. deviation 6fit < 0.7%, while the LLD fits can leave a deviation 
of as much as 3% (in the case of the B profiles). The LLD fits are marginally better than 
the PCA fits only for some of the R3750 and R4000 profiles. The more striking difference 
is in the excess flux. While the fluxes of the PCA flts differ from the actual source fluxes 
by less than 2 x 10~^, the LLD fits overestimate the source flux for all B, V, R3500 and 
R3750 profiles, by as much as 0.75% . This renders the LLD model also less useful for fitting 
microlensing light curves, where such a difference is not negligible. In the remaining profiles 
the flux is underestimated; a flux accuracy comparable to the PCA flts is never achieved. 

Although our PCA basis was constructed from limb darkening proflles of cool giant 
models. Figure 11 demonstrates this PCA model can be better suited than the LLD model 
even for fltting measured data of narrow-band solar limb darkening. For this flgure we 
used data published by Mitchell (1959) spanning the BVRI spectral range, namely his 
measurements at 450, 545, 660 and 848 nm. The radial coordinates of the data points are 
marked in the figure by crosses. In order to give even weight to all radii^^, we interpolated 
the observed data by a spline in the manner described at the beginning of this appendix 
to obtain B{r). We performed least-squares flts to these interpolations using both limb 
darkening models and plotted the deviation curves B{r)/Bfit{r) — 1. The small- amplitude 
bumps in the curves are due to the interpolation. As seen from the flgure, while for 450 nm 



'Computed as unweighted fits to the spHne-interpolated profiles (see Fig. 2). 

Note that the narrow region close to the limb has the highest density of data points. 
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and 848 nm both models have a comparable maximum deviation (within the range r < 0.9), 
for 545 nm and 660 nm the PCA model fits the data better. It is interesting to note that 
here the LLD model systematically provides fits more peaked than the original profiles. 

A more rigorous approach for studying different stellar types would be to construct 
PCA bases from model atmospheres of such stars. The reason for the cool giant PCA basis 
fitting the solar data is primarily physical. For both stellar types as well as for a wide range 
of other cool stars, the dominant source of opacity determining the limb darkening in the 
optical range is the same — namely, the H~ ion (e.g.. Gray 1992). 



-33- 



REFERENCES 

Afonso, C, et al. 2000, ApJ, 532, 340 
Albrow, M. D., et al. 1999a, ApJ, 522, 1011 
Albrow, M. D., et al. 1999b, ApJ, 522, 1022 
Albrow, M. D., et al. 2000, ApJ, 534, 894 
Albrow, M. D., et al. 2001, ApJ, 549, 759 
Alcock, C, et al. 1997, ApJ, 491, 436 
An, J. H., et al. 2002, ApJ, 572, 521 
Becker, A., et al. 1997, BAAS, 29, 1347 
Bessell, M. S. 1990, PASP, 102, 1181 

Bogdanov, M. B., & Cherepashchuk, A. M. 1995, Astron. Rep., 39, 779 

Bogdanov, M. B., & Cherepashchuk, A. M. 1996, Astron. Rep., 40, 713 

Bogdanov, M. B., & Cherepashchuk, A. M. 1999, Astron. Rep., 43, 601 

Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, R. W., & Burrows, A. 2001, ApJ, 
552, 699 

Bryce, H. M., & Hendry, M. A. 2000, to appear in ASP Conf. Ser., Microlensing 2000, ed. 
J. W. Menzies & P. D. Sackett (San Francisco: ASP), preprint (astro-ph/0004250) 

Burns, D., et al. 1997, MNRAS, 290, Lll 

Byrd, P. F., & Friedman, M. D. 1971, Handbook of Elliptic Integrals for Engineers and 
Scientists (2nd ed.; Berlin: Springer) 

Craig, I. J. D., & Brown, J. C. 1986, Inverse Problems in Astronomy (Bristol: Adam Hilger) 

Cummings, I. N., Hearnshaw, J. B., Kilmartin, P. M., & Gilmore, A. C. 1998, in ASP Conf. 
Ser. 135, A Half-Century of Stellar Pulsation Interpretations, ed. P. A. Bradley & J. 
A. Guzik (San Francisco: ASP), 213 

Deeg, H. J., Garrido, R., & Claret, A. 2001, New Astronomy, 6, 51 

Gaudi, B. S., & Gould, A. 1999, ApJ, 513, 619 



-34- 



Gaudi, B. S., & Sackett, P. D. 2000, ApJ, 528, 56 

Gould, A. 2000, ApJ, 542, 785 

Gould, A., & Welch, D. L. 1996, ApJ, 464, 212 

Gray, D. F. 1992, The Observation and Analysis of Stellar Photospheres (2nd ed.; Cambridge: 
Cambridge University Press) 

Gray, N., & Coleman, 1. J. 2000, to appear in ASP Conf. Ser., Microlensing 2000, ed. J. W. 
Menzies & P. D. Sackett (San Francisco: ASP), preprint (astro-ph/0004200) 

Griest, K., & Hu, W. 1992, ApJ, 397, 362 

Griest, K., & Safizadeh, N. 1998, ApJ, 500, 37 

Grygar, J., Cooper, M. L., & Jurkevich, I. 1972, Bull. Astr. Inst. Cz., 23, 147 
Hajian, A. R., et al. 1998, ApJ, 496, 484 
Han, C, & Gould, A. 1997, ApJ, 480, 196 
Han, C, & Jeong, Y. 1998, MNRAS, 301, 231 

Hendry, M. A., Coleman, I. J., Gray, N., Newsam, A. M., & Simmons, J. F. L. 1998, NewAR, 
42, 125 

Heyrovsky, D., & Loeb, A. 1997, ApJ, 490, 38 

Heyrovsky, D., & Sasselov, D. 2000, ApJ, 529, 69 

Heyrovsky, D., Sasselov, D., & Loeb, A. 2000, ApJ, 543, 406 (HSL) 

Jacob, A. P., Bedding, T. R., Robertson, J. G., & Scholz, M. 2000, MNRAS, 312, 733 

Jha, S., Charbonneau, D., Garnavich, P. M., Sullivan, D. J., Sullivan, T., Brown, T. M., & 
Tonry, J. L. 2000, ApJ, 540, L45 

JoUiffe, I. T. 1986, Principal Component Analysis (New York: Springer) 

Jorissen, A., Mowlavi, N., Sterken, C, & Manfroid, J. 1997, A&A, 324, 578 

Koen, C, & Laney, D. 2000, MNRAS, 311, 636 

Loeb, A., & Sasselov, D. 1995, ApJ, 449, L33 



-35- 



Mitchell, W. E. 1959, ApJ, 129, 93 
Mozurkewich, D., et al. 1991, AJ, 101, 2207 
Paczynski, B. 1996, ARA&A, 34, 419 
Percy, J. R., & Parkes, M. 1998, PASP, 110, 1431 

Percy, J. R., & Polano, S. 1998, in ASP Conf. Ser. 135, A Half-Century of Stellar Pulsation 
Interpretations, ed. P. A. Bradley & J. A. Guzik (San Francisco: ASP), 249 

Popper, D. M. 1984, AJ, 89, 132 

Popper, D. M., Andersen, J., Clausen, J. V., & Nordstrom, B. 1985, AJ, 90, 1324 

Preisendorfer, R. W. 1988, Principal Component Analysis in Meteorology and Oceanography 
(Amsterdam: Elsevier) 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes 
in C (2nd ed.; Cambridge: Cambridge University Press) 

Quirrenbach, A., Mozurkewich, D., Buscher, D. F., Hummel, C. A., & Armstrong, J. T. 
1996, A&A, 312, 160 

Rhie, S. H., et al. 1999, ApJ, 522, 1037 

Rhie, S. H., & Bennett, D. P. 1999, ApJ, submitted, preprint (astro-ph/9912050) 

Sasselov, D. D. 1996, in ASP Conf. Ser. 109, Cool Stars 9, ed. R. Pallavicini & A. K. Dupree 
(San Francisco: ASP), 541 

Simmons, J. F. L., Newsam, A. M., & Willis, J. P. 1995, MNRAS, 276, 182 

Uitenbroek, H., Dupree, A. K., & Gilliland, R. L. 1998, AJ, 116, 2501 

Valls-Gabaud, D. 1995, in Large Scale Structure in the Universe, ed. J. P. Miicket, S. 
Gottlober, & V. Miiller (Singapore: World Scientific), 326 

Witt, H. J. 1995, ApJ, 449, 42 



This preprint was prepared with the AAS macros v5.0. 



-36- 




Fig. 1. — Dashed line: Angle-integrated amplification A{l,r) for a lens with Einstein radius 
e=l as a function of radial distance r for lens positions /=0, 0.25, 0.5, 0.75 and 1 (curves with 
peaks from left to right). Solid line: A{l,r) / e for a lens with e=10, same lens positions. 




Fig. 2. — Normalized BVRI brightness profiles of three logg — 0.5 model atmospheres (from 
upper to lower at r = in each filter): T=3500, 3750, 4000 K. Radial positions of computed 
points are marked at bottom, plotted curves are spline interpolations to these points. 



-38- 




Fig. 3. — First three orthonormal functions of the PCA basis for brightness profile decompo- 
sition, obtained by principal component analysis of 32 brightness profiles of red giant model 
atmospheres (see Table 1). 




Fig. 4. — Recovered PCA (upper panels) and LLD (lower panels) limb darkening parameters 
for model profiles: B3750gl0 (left column), R3500g05 (central column) and I4000gl0 (right 
column). Individual points correspond to values of the original Einstein radius eo and impact 
parameter po directly above them in the table, triplets of neighboring points correspond to 
blend parameters Po = 0, 0.25, 0.5 . Dotted vertical bars: range of values from all simulations. 
Solid horizontal lines: original best-fit values. Dotted horizontal lines: full parameter range. 
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Fig. 5. — Brightness profile inversion accuracy from center to limb of lensed star, for cq = 10 . 
Bold solid line: bias 5b \ thin solid line: "median median" 5mm \ bold dashed lines: total 
errorbars; thin dashed lines: systematic errorbars. Dotted lines: 5 = 0, ±2% and r = 0.9. 
Multiply total errorbars by inverted profile to obtain absolute errorbars. Columns correspond 
to marked impact parameters po- Upper panels: PCA inversion; lower panels: LLD inversion. 
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Fig. 6. — Baseline detail of the MACHO B (upper panel) and MACHO R (lower panel) light 
curves of the MACHO Alert 95-30 event. Vertical axes correspond to relative flux devia- 
tions from the best-fit microlensing solution baselines (marked by dotted lines). Note the 
correlated variability as well as the systematically positive MACHO B deviation. Following 
Alcock et al. (1997), time is measured from JD 2,448,623.50 . 
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Fig. 7. — Relative residuals Fi/ F{li) — 1 of the best-fit microlensing solution close to the 
peak of the MACHO Alert 95-30 event. Individual panels correspond to hght curves marked 
in the left corners. The time of closest approach and limb crossing times of the best-fit 
solution are marked by vertical dotted hues (time is measured from JD 2,448,623.50) . 
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Fig. 8. — Lower panel: BVRI amplification light curves (from top to bottom at t = 0) of 
a microlensing transit of a T=3750 K, logg = 0.5 red giant by an e = 10 lens with a zero 
impact parameter. Upper panel: corresponding brightness profiles normalized to unit flux 
(same sequence at t — 0), horizontal scale matches the lens position. 
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Fig. 9. — Broadband microlensing chromaticity: relative variation of microlensing amplifi- 
cation from Figure 8 with filter passband, as a function of lens position (in source radius 
units). Curves correspond to Einstein radn e = 0.2, 0.5, 1, 2, 5, 10 (from lowest curve). 
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Fig. 10. — Full parameter range of limb darkening profiles using the 2-term PCA (left panel) 
and LLD bases (right panel), normalized to unit flux and plotted on same scale as Figure 2. 
Limb darkening parameter values k are spaced by 0.03 in the PCA model, by 0.1 in the 
LLD model. Solid curves correspond to /t = in each model (first basis term). The vertical 
dotted line marks the common intersection point, characteristic of each model. 
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Fig. 11. — Comparison of PCA (solid line) and LLD (dashed line) residuals of fits to measured 
solar limb darkening data (Mitchell 1959) at four selected wavelengths as denoted. Radial 
positions of the data points are marked at the lower edge of the plots. Dotted hnes serve for 
orientation: B/Bfn — 1 = 0, ±1% and r = 0.9 . 
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Table 1: List of model brightness profiles used in this work. Adopted profile names indicate 
the filter, effective temperature and surface gravity of the model atmosphere. Following 
columns: best-fit PCA and LLD limb darkening parameters k, normalized r.m.s. deviations 
Sfit and relative excess fluxes AF/F of the flts. 



